Three-dimensional MgB2-type superconductivity in hole-doped diamond. 
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We substantiate by numerical and analytical calculations that the recently discovered superconductivity below 
4 K in 3% boron-doped diamond is caused by electron-phonon coupling of the same type as in MgB2, albeit 
in 3 dimensions. Holes at the top of the zone-centered, degenerate a-bonding valence band couple strongly to 
the optical bond-stretching modes. The increase from 2 to 3 dimensions reduces the mode-softening crucial for 
Tc reaching 40 K in MgB2. Even if diamond had the same bare coupling constant as MgB2, which could be 
achieved with 10% doping, Tc would only be 25 K. Superconductivity above 1 K in Si (Ge) requires hole-doping 
beyond 5% (10%). 

PACS numbers: 74.70.-b,74.70.Ad,74.25.Kc 



Recently, superconductivity below Tc ~4 K was reported in 
diamond doped with x ^3% boron, that is, with ^0.03 holes 
per carbon atom [IJ. Such high hole-doping levels can be 
achieved due to the small size of boron. It had previously been 
observed that the prominent Raman line caused by the zone- 
center optical phonons at 1332 cm^^ downshifts and broadens 
significantly upon heavy boron doping j^l • In this Letter we 
shall make plausible that the superconductivity in hole-doped 
diamond is due to the coupling of the holes to the optical zone- 
center phonons, a mechanism similar to the one causing high- 
temperature superconductivity in MgB2, but without some of 
its interesting features. We shall also estimate transition tem- 
peratures for hole-doped Si and Ge. 

The discovery f?) of superconductivity below 40 K in 
MgB2, a binary compound isostructural and isoelectronic 
with graphite, came as a surprise for the scientific community. 
By now, it is well understood what the mechanism is and why 
MgB2 is special |4, 5, 6, 7, 8|: In contrast to other known 
sp^ -bonded superconductors, such as intercalated graphite, 
alkali doped fullerides, and organic superconductors whose 
charge carriers are exclusively 7r-electrons, MgB2 has holes 
at the top of the bonding cr-bands at the zone center. These 
holes, on two narrow Fermi cylinders with radii ^--^1/5 of the 
Brillouin-zone radius, couple strongly (A ^ 1) to the two op- 
tical bond-stretching modes with q < 2fci? <C kszi giving 
rise to a strong 2-dimensional Kohn anomaly in the phonon 
spectrum. This strong coupling between a few zone-center 
holes and optical phonons is what drives the high-temperature 
superconductivity in MgB2. Experience shows |5|, and it can 
be proved for parabolic bands with 2fci? ^ ksz |9], that the 
coupUng constant is given by the Hopfield expression. 



where N is the density of states (DOS) per spin at the Fermi 
level of the a holes. Moreover, ±_Du is the splitting of the de- 
generate top of the (T-band by the displacement eu of a frozen, 
optical zone-center phonon with normalized eigenvector e and 
energy lo. The optical phonons are softened by the interac- 



tion with the holes, cj^ ~ cjq/ (1 + 2A) when q < 2kp, and 
that significantly enhances A and Tc ^ cj cxp (— 1/A). This 
softening is presumably weakened by anharmonicity I loi fllll . 
The DOS is independent of doping because the cr-band is 2- 
dimensional. As a consequence, a decrease in the number of 
holes, e.g. by carbon-doping, should not cause A to decrease, 
except through the anharmonic hardening of u caused by the 
decrease of Ep [12]. In stoichiometric MgB2 there are more 
carriers in the tt bands than in the a bands (0.09 per B), but 
the former couple far less to phonons than the latter, and since 
there seems to be very little impurity scattering between the a- 
and TT-bands, MgB2 is the first superconductor which clearly 
exhibits multiple gaps below a common Tc 10, S [13. Ell • 

Instead of having 7r-bands and three 2-dimensional bond- 
ing cr-bands, sp'^-bonded semiconductors like diamond have 
four 3-dimensional bonding cr-bands. The top of this valence 
band is three-fold degenerate with symmetry T2g, and so are 
the zone-center optical phonon modes. Like in MgB2 , cr-holes 
with small k p should therefore couple strongly, and for small 
kp exclusively, to the optical bond-stretching modes, with the 
main differences being that in 3 dimensions the Kohn anomaly 
is weaker and the DOS increases with hole doping like kp, 
the radius of the average Fermi sphere. Since there are 3 
bands and two carbon atoms per cell, (kp/ksz)'^ — x/3. For 
x=0.03, kp/ksz is 0.22, which is like in MgB2. Due to the 
lack of a metallic 7r-band, diamond becomes an insulator once 
X is below 1-2 per cent. We shall now substantiate this sce- 
nario for the observed superconductivity in hole-doped dia- 
mond by providing quantitative details, and we shall also con- 
sider the possibility of superconductivity in hole-doped Si and 
Ge. In particular, we shall present results of density-functional 
(LDA) calculations and estimate Tc using Eliashberg theory. 

A substitutional boron impurity in diamond has an accep- 
tor level with binding energy 0.37 eV |14]. With increased 
doping, the boron impurity band will eventually overlap the 
diamond valence band and the system becomes metallic at a 
boron concentration of 8 • 10^'^ cm^'^. Since this is one order 
of magnitude lower than the doping at which superconductiv- 
ity was observed 1 1], we felt justified in using a virtual crystal 
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approximation (VCA) in which the carbon nuclei have charge 
Q — X and the crystal is neutral. 

The valence bands were calculated with the scalar- 
relativistic full-potential LMTO method fl5], and the phonon 
dispersions and the electron-phonon spectral function a^F 
were calculated with the linear-response method 1 15|. Effects 
of anharmonicity were considered in a second step. We used a 
triple-K spd LMTO basis set and represented the charge den- 
sities and potentials by spherical harmonics with 2 < 6 in- 
side non-overlapping muffin-tin spheres and by plane waves 
with energies less than 400 Ry between the spheres. The re- 
sulting band structure for undoped diamond agrees with those 
of earlier LDA calculations. Due to the smallness of kp, we 
needed to use a fine k-mesh chosen as a 1/32^ sublattice in 
reciprocal space. a^F is evaluated as a weighted sum over 
linewidths of individual phonons, and for this a fine, yet af- 
fordable q-mesh is needed. It was chosen as a 1 sublattice 
in reciprocal space. The A-values obtained herewith are ac- 
curate when X > 0.05, whereas Eq. Q which we have now 
derived analytically also in three dimensions, is more accurate 
for smaller dopings. Since, even for 10% doping, we calcu- 
lated an increase of the lattice constant by less than a per cent, 
we did use the experimental lattice parameters for the Mndoped 
materials in all subsequent calculations . 

In Fig.^we show the top of the valence-band structure cal- 
culated for 10% hole-doped diamond. For this unrealistically 
heavy doping, N reaches 75% of the a-band DOS in MgB2. 
The electronic parameters may be found together with those 
for MgB2 in the first columns of the table. Due to the devia- 
tions from parabolicity seen in Fig.ffl N decreases somewhat 
faster than x^^"^. As is well known, the LDA gap is too small 
and this leads to a slight underestimate of the valence-band 
masses and the DOS. Nevertheless, properties derived from 
the total energy, like phonon energies, are quite accurate. 

For the displacement eu of a frozen, optical zone-center 
phonon with two bonds stretched and two contracted, say 
those in respectively the positive and negative z-directions, 
the top of the valence band is deformed as shown in Fig. ^ 
For small k there are two identical bands, split in energy by 



±Du, and a band which does not move with respect to the 
Fermi level. At F the corresponding wavefunctions are those 
linear combinations of the bond orbitals which have respec- 
tively (px ± py ) \/2 and p^ symmetry. The value of D given 
in the table is seen to be nearly twice as large as in MgB2. For 
pure diamond it agrees with the accepted value lUiil . and it is 
seen to decrease slightly with doping. 

In Fig. 12] we show the phonon dispersions calculated in 
the harmonic approximation for undoped and hole-doped di- 
amond. The dispersions for pure diamond, including the 
slight upturn of the uppermost mode when moving away from 
the zone center 1 17), are well reproduced, and for the fre- 
quency 1332 cm^^ of the optical zone-center modes we calcu- 
late ti)o=1292 cm^^. Previous LDA calculations 1 18] obtained 
similar results. In the presence of hole-doping, the calculated 
dispersions of the optical modes clearly exhibit softening near 
the zone center and a 3 -dimensional Kohn anomaly around 
q — 2kF- 

The softening of the zone-center phonons is (2/3) A instead 
of A, as in the case of MgB2. This is most easily seen by con- 
sidering a frozen phonon calculation and Fig.Q The adiabatic 
redistribution of (1/3) 2NDu electrons from the upper third 
to the lower third of the deformed valence band decreases the 
energy of each electron by Du, and therefore perturbs the po- 
tential energy of the harmonic oscillator by — (1 /3) 2ND^v?. 
As a consequence, {1/2) Muj'^ = (1/2) Afcj^ - (1/3) 2iVD2, 
and by use of Eq. Q we get: uj"^ = t^o/(l + 2|A). In MgB2 
no part of the cr-band is passive in the screening of the phonon, 
so the factor 2/3 is missing. The value of A deduced from the 
frequencies, u) and tjQ = io{x—Q), of the optical zone-center 
modes calculated by linear response is given in the table (A^^). 
This A[^-value is seen to agree well with the value Ad obtained 
by use of Eq. ([Q. In order to separate the materials and di- 
mensional dependencies of A, we express it in terms of a bare 
coupling constant, Aq , and an enhancement due to the phonon 
softening: 
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FIG. 1: LDA band- structure of diamond with a;=0.1 holes/C (full 
lines). A frozen optical zone-center phonon with two bonds stretched 
and two contracted by u/\/2 splits the bands by ±Du (dotted lines). 
u=0.05A. 



The enhancement is weaker in 3 dimensions than in 2, where 
the reduction factor 2/3 is missing. As for the materials depen- 
dence, the Ao-values given in the table first of all show that 
10% doped diamond has the same Ao ^ 1/3 as MgB2: The 
bare force constant, Mwg, is 0.49 times its value in MgB2, N 
is 0.75, and D is 1.65. Due to the difference in dimensional- 



ity, A 



1/3 



~ 0.6 in doped diamond, but A ^ 1-2/^ — 1 
in MgB2. With decreasing doping in diamond, N decreases 
roughly like x^/'^, D increases slightly, and Muj"^ is constant. 
As a consequence, for 3% doping Ao is only 0.21 and A is 
0.30. 

We can also calculate the electron-phonon spectral func- 
tion and A = 2 ^ uj^^ F {uj) dio numerically hy sampling 
over all phonon branches and energy bands. The result shown 
in Fig. |3] confirms that only the optical phonons interact with 
the holes: a^F vanishes for phonon frequencies below that 
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TABLE I: N is in states/eV/spin/f.u. D is in eV/A. uj is in cm~ ^ . Aq is tiie bare electron-piionon coupling constant defined in Eq. J2j- -^d 'T"! 
Ao; are estimates of the coupling constant as obtained from respectively Eq. and the softening of lo. A is obtained from the numerical linear- 
response calculation and includes all phonons and tj-electrons; for MgB2 it is Xca Ii21.1 . a = E'p/Dt. \a is A corrected for anharmonicity 
using Eq.js}. Tc is obtained from Eq. J4j using Aa, cOa, and ^*=0.1. 
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of the optical zone-center modes, then jumps to a maximum, 
and finally falls. The decay occurs more slowly than in MgB2 
due to the increase of dimensionality. The A-values (A) ob- 
tained from this calculation again agree well with those ob- 
tained from Eq. Q and from the phonon softening. 

In MgB2 the role of anharmonicity of the optical phonon 
modes with q < 2kF has been stressed |10]; it hardens the 
phonon by about 20% and thus decreases A from 1.0 to 
Aa=0.78, as given m the table. However, this has recently 
been questioned |11]. While anharmonicity may be crucial 
in MgB2 it has at most a noticeable effect on the supercon- 
ductivity in diamond at small dopings (x <3%) as we shall 
now see: The anharmonicity appears in frozen phonon cal- 
culations (see Fig. 0, because once the displacement u ex- 
ceeds Ep/D, the lower band is full so that the screening is 
lost 1 12]. In the expression for the perturbation of the poten- 
tial energy of the oscillator, should now be substituted by 
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FIG. 2: Phonon dispersions calculated with the linear-response 
method for diamond with a; = 0, 1, 5, and 10 % hole doping. 



{\u\ - Ep/D) e [Ep/D - \u\), provided that we simplify 
the DOS shape by a square. It has been shown that the most 
important anharmonic contribution to Tc is the decrease of the 
first excitation energy 1 19]. By first-order perturbation the- 
ory, th is is sim ply (1/3) 2NDH'^ [1 - erf (Ep/Dt)] , where 
UJ is the classical turning point in the ground state. 
Introducing again Eq.([l) we obtain the result: 
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A ^ tj2 - 1 + 2 (2/3) A [1 - erf {Ep/Dt)] ' 



(3) 



For MgB2 the assumption of a square N (E) is good, but due 
to the missing factor 2/3 and the presence of the 7r-band, (2/3) 
in Eq. Q should be substituted by [1 - N^/ {N + 2N^)]. 
For hole-doped diamond, N{E) has square-root shape, and 
this we crudely take into account by substituting Ep in Eq. 
^ hy E'p = (2/3) Ep. In the table we have included the 
ratio Ep /Dt = a as well as the results for Aa. We see that the 
effect of anharmonicity may be important in MgB2 but merely 
noticeable in hole-doped diamond Eoil . 

For this type of superconductivity, which is characterized 
by an Eliashberg function with the shape exhibited in Fig.|3l 
and which we can idealize by a 5-function at the frequency oj 
of the optical zone-boundary phonon, solution of the Eliash- 
berg equations yields with high accuracy: 



Lu exp 



(4) 



This is McMillan's expression with all numerical factors, 
which he obtained by fitting to F (uj) of niobium, set equal 
to unity. For the cases considered in the present paper, it does 
not make much difference whether one uses McMillan's fac- 
tors or unity inside the exponential, but it is important that the 
prefactoris uj, rather than {u!\n) /1. 2. 

We can finally estimate Tc from Eq. Q with the values for 
Aa given in the table and uja from Eq. (|3}- For the Coulomb 
pseudopotential, the standard value /i*=0.1 was used in all 
cases. For MgB2 we neglected the 7r-bands. Considering the 
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pie analytical expressions for MgB2-type superconductivity. 

A purely electronic mechanism for the observed supercon- 
ductivity was recently suggested |22|, and after submission 
of the present manuscript two works similar to ours appeared 
1E31 13 . The latter used a supercell approach to simulate the 
boron doping and found an electron-phonon coupling in very 
good agreement with our results. 

We are grateful to O. Dolgov, M. Cardona, G. B. Bachelet, 
E. Cappelluti, and L. Pietronero for many interesting discus- 
sions. 



FIG. 3: The phonon density of states, F{uj) , and Eliashberg function, 
a^Fliu), calculated numerically by linear response considering all 
electrons and phonons. 

uncertainties in our calculation of A and uj, the uncertainty of 
/i*, and the experimental estimation of the doping level, we 
do find critical temperatures in good agreement with present 
experimental knowledge. We therefore believe to have sub- 
stantiated our claim that the superconductivity in hole-doped 
diamond is of MgB2-type, but in three dimensions. 

We repeated our calculations for hole-doped Si and Ge, and 
include those results in the table for which Ep largely ex- 
ceeds the spin-orbit splitting, which we neglected. Whereas 
hole-doped C shows superconductivity above 1 K for doping 
levels presently obtainable. Si and Ge seem to need twice as 
high doping levels. The main reason is that the deformation 
potential in Si and Ge is about four times smaller than in C, 
which is too small to take advantage of having twice as large 
a DOS and a three times smaller force-constant. There is also 
a qualitative difference to diamond: For heavily doped Si and 
Ge, the holes not only couple to the optical, but also to the 
acoustic phonons. This is the reason why A exceeds Xd ^ A^^ . 

In conclusion, we have shown that the recently discov- 
ered superconductivity in hole-doped diamond below 4K is 
of MgB2-type, but in three dimensions. This means that 
the mechanism is coupling of a few holes at top of the cr- 
bonding valence band to the optical bond-stretching zone- 
center phonons. The increase from 2 to 3 dimensions limits 
the strong softening of the optical modes mainly responsible 
for the high Tc in MgB2. On the other hand, the deforma- 
tion potentials in diamond are twice stronger than in MgB2. 
Kelvin-range superconductivity in Si and Ge would require 
hole-doping levels of 5-10%. Finally, we have obtained sim- 
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